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Some important aspects of the electronic structure of the iron oxypnictides depend very sensitively 
on small changes in interatomic distances and bond angles within the iron-pnictogen subunit. Using 
first-principles full-potential electronic structure calculations, we investigate this sensitive depen- 
dence, contrasting in particular LaOFeAs and LaOFeP. The width of the Fe-bands is significantly 
larger for LaOFeP, indicating a better metal and weaker electronic correlations. When calculated 
at their experimental crystal structure these two materials have significantly different low-energy 
band structure. The topology of the Fermi surface changes when going from LaOFeP to LaOFeAs, 
with a three-dimensional hole pocket present in the former case transforming into a tube with two- 
dimensional dispersion. We show that the low-energy band structure of LaOFeAs evolves towards 
that of LaOFeP as the As atom is lowered closer to the Fe plane with respect to its experimental 
position. The physical origin of this sensitivity to the iron-pnictogen distance is the covalency of the 
iron-pnictogen bond, leading to strong hybridization effects. To illustrate this, we construct Wannier 
functions, which are found to have a large spatial extension when the energy window is restricted to 
the bands with dominant iron character. Finally, we show that the Fe bandwidth slightly increases 
as one moves along the rare-earth series in ReOFeAs and discuss the physical origin of this effect. 

PACS numbers: 



Introduction 

The recent discovery of superconductivity in electron- 
doped rare-earth oxyarsenidesA has generated a great 
deal of interest in the electronic structure and magnetic 
properties of these compounds, with different types of 
magnetic 2 ^ 3 - or charge 4 - fluctuations conjectured to be at 
the origin of the pairing mechanism. 

The rare-earth iron oxyarsenides belong to a wider 
class of rare-earth oxypnictides ReOTPn (where Re and 
T are the rare-earth and transition metals, respectively, 
Pn = P, As is the pnictogen ion, with nominal charge —3) 
which have the same tetragonal structure of ZrCuSiAs- 
type&. In spite of obvious similarities in the overall elec- 
tronic structure, those compounds differ in their mag- 
netic properties and their tendency towards supercon- 
ductivity. It is particularly interesting in this respect to 
contrast the properties of LaOFeAs to those of its phos- 
phate homologue LaFeOP. Relating the observed differ- 
ences to relevant differences in their electronic structure 
may provide important clues into the nature and the ori- 
gin of magnetism and superconductivity in these mate- 
rials. The undoped stoechiometric compounds have dif- 
ferent magnetic properties: LaOFeP is non-magnetio^, 
while LaOFeAs undergoes a magnetic ordering transi- 
tion at a temperature of 134 K— preceded by a struc- 
tural transition at a slightly higher temperature^. Their 
transport properties above the magnetic transition tem- 
peratures are also different, with the electrical conduc- 
tivity of LaOFeP being substantially larger than that 
of LaOFeAs 1 -^. These compounds also differ signifi- 



cantly in their superconducting properties. While un- 
doped LaOFeAs is magnetic and non-superconducting, 
the situation in the absence of doping for the parent com- 
pound LaOFeP is currently somewhat controversial. A 
superconducting transition in the range 4-7 K was orig- 
inally reported^, and recently confirmed in Ref. 0, but 
questioned in Ref. [1(3 which found no superconductivity 
above 0.35 K. At any rate, electron doping suppresses the 
antiferromagnetism in LaOFeAs and leads to a supercon- 
ducting phase with T c as high as 26 K, while it does not 
induce such drastic changes in the phosphate compound, 
in which it causes only a moderate increase of T c up to 
8 K&. 

The electronic structure of the LaOFeP compound, 
in particular its Fermi surface (FS), has been investi- 
gated theoretically by LebegueA 1 -. His results can be 
outlined as follows: the LaOFeP FS consists of two 
ellipsoidal electronic cylinders centered at the M — A 
line at the Brillouin zone (BZ) edge, and of two hole 
cylinders centered at the BZ center T-Z line, together 
with a single hole pocket with three-dimensional (3D)- 
like dispersion at the Z point. A basically identical 
shape of the FS has been proposed in several theoretical 
paper a 12 i 13 i 14 for LaOFeAs , with electron doping result- 
ing in the disappearance of the 3D-hole pocket, which 
gets filled, thus changing the FS topology and increasing 
the two-dimensional character upon doping. This topo- 
logical change and reduction of effective dimensionality 
has been suggested to be at the origin of the different 
types of magnetic fluctuations observed in the stoichio- 
metric and doped LaOFeAs compounds, respectively 3 ^ 4 -, 
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and possibly also responsible for the appearance of su- 
perconductivity in the doped compound. It is impor- 
tant to note, in this respect, that most electronic struc- 
ture calculations in Refs i 12 i 13 i 14 have been carried out 
at values of the lattice constants and internal parameters 
(the Fe-As and La-0 interplane distances) which were ob- 
tained by relaxing the structure within LDA/GGA. As 
pointed out very recently by Mazin et al/H, rather small 
changes in the LaOFeAs structural parameters, partic- 
ularly, in the Fe-As interatomic distance, may lead to 
the disappearance of the 3D hole pocket from the FS of 
stoichometric LaOFeAs. Furthermore, internal parame- 
ters of the LaOFeAs experimental crystal structure are 
rather poorly reproduced by theoretical LDA/GGA cal- 
culations, with the Fe-As interplane distance being un- 
derestimated by about 6%^. 

For these reasons, a comparative study of the 
LaOFeP and LaOFeAs compounds is especially relevant 
in order to clarify the relations between electronic prop- 
erties and the electronic structure of these materials. We 
do know that stoichometric LaOFeP and LaOFeAs differ 
substantially, particularly in properties directly related 
to the low-energy electronic structure such as magnetic 
ordering, transport or superconductivity. Therefore, one 
may hope and expect that electronic structure calcula- 
tions carried out at judiciously chosen structural param- 
eters should be able to capture this difference. In this ar- 
ticle, we carry out a detailed first-principles investigation 
of LaOFeAs and LaOFeP aimed at identifying relevant 
differences in their electronic structures while relating 
them to the small but significant differences which ex- 
ist between the LaOFeAs and LaOFeP lattice structures. 
We show that electronic structure calculations at the ex- 
perimental lattice parameters and internal positions re- 
sult in the Fermi surfaces of LaOFeP and LaOFeAs hav- 
ing different shapes, with the 3D hole pocket present in 
the case of LaOFeP , but not in the case of LaOFeAs. 
These changes in the FS topology are due to shortening 
of the iron-pnictogen bond-length in LaOFeP as com- 
pared to LaOFeAs. We identify the sensitivity of the elec- 
tronic structure to this bond-length as being due to the 
rather high degree of covalency associated with the iron- 
pnictogen bond, and illustrate this point by constructing 
appropriate Wannier functions. 

The paper is organized as follows. We briefly outline 
our calculational approach in Sec. |H The overall elec- 
tronic structure of LaOFeAs and LaOFeP is discussed in 
Sec. [ill Low-energy aspects of this electronic structure 
and Fermi surfaces are described in Sec. IIIIl and the 
differences between the two materials are illustrated by 
a study of the sensitivity of these low-energy aspects to 
the Fe-pnictogen distance. Then, in Sec. IIV1 we per- 
form a Wannier function construction which illustrates 
this sensitivity and relates it to the covalency of the iron- 
pnictogen bond. Finally, in Sec. |V]we study the influence 
of the structural properties on the width of the Fe- bands 
and its material dependence. 



Table I: Crystal-structure data for LaOFeAs^ and LaOFeP 7 . 
Lattice constants a and c, internal positions of the Pn and 
La (z) in units of c, and the vertical distance of the Pn with 
respect to the Fe plane (d v ). The vertical distance is related 
to the z-parameter by d v = (z— 0.5)c, and the angle 9 between 
the As-Fe bond and the Fe-plane by tan(8) — (2z — l)c/a. 





LaOFeAs 


LaOFeP 


a 


4.0301 A 


3.9636 A 


c 


8.7368 A 


8.5122 A 


z (La) 


0.1418 


0.1487 


z{Pn) 


0.6507 


0.6339 


d v 


1.32 A 


1.14 A 


9 


33.2° 


29.9° 


dFe-As 


2.41 A 


2.33 A 



I. COMPUTATIONAL DETAILS 

The first principle calculations of the band-structure 
and Fermi surface of these materials presented in Sec. 
II and Sec. Ill were performed using the Full-Potential 
APW+local orbitals method as implemented in the 
Wien2k 1 - code. We considered 800 k-points in the BZ 
and checked that all the properties presented in this ar- 
ticle were converged with this mesh. The results that we 
present here have been obtained within the Local Den- 
sity Approximation (LDA)i£ to the exchange-correlation 
potential but it is worth mentioning that our conclusions 
do not vary if we use the generalized-gradient approxima- 
tion (GGA)— instead. We have also performed a calcula- 
tion using the Quantum-ESPRESSO package^, with the 
computational details like in [2(J and constructed maxi- 
mally localized Wannier functions^ by using the code 
developed by Mostofi et al.— . Last but not least, we have 
performed an Nth order muffin tin orbital construction 
for the Fe-d orbitals^. 



II. BAND STRUCTURE OF LaOFeAs AND 
LaOFeP : COMPARISON ON A LARGE 
ENERGY-SCALE 

In this section, we compare the calculated band struc- 
ture of LaOFeP and LaOFeAs at their experimental crys- 
tal structure, on a large energy scale covering the whole 
Fe, Pn (=As or P) and O bands as shown in Fig. [TJ 
The crystal structure parameters were taken from exper- 
iments and are given in Table [H 

The obtained bandstructures for both materials are de- 
picted on Fig.[TJ In this figure, the colors of the bands are 
simply guides to the eyes in order to compare more eas- 
ily their relative positions and bandwidths. The bands 
plotted in red and green have mainly Fe 3d and Pn-0 p 
character, respectively. The bands shown in black corre- 
spond mainly to La 4/ states which are well above the Fe 
bands. The blue arrows in between the two plots show 
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the approximate bandwidth of the Fe 3d bands in each 
compound. One may observe that the bands that are 
close to the Fermi level, which correspond to Fe, have 
larger bandwidth in LaOFeP as compared to LaOFeAs. 
This is a plausible explanation of why the former is ob- 
served to be a better metal (in the sense of smaller resis- 
tivity This also indicates a lesser degree of electronic 
correlations in LaOFeP. 

The difference in bandwidth results from a combina- 
tion of three factors. First, from the difference in Fe-Fe 
distance which is significantly smaller in LaOFeP. Sec- 
ondly, P and As have different electronegativity, which 
leads to a shift towards more negative energies of the P- 
bands in LaOFeP. Finally, the Pn has different vertical 
distances with respect to the Fe plane, a quantity which 
is mainly governed by the internal z parameter, and may 
influence the width of the Fe' s bands through the differ- 
ent strength of the Fe-Pn bond. 



LaOFeAs LaOFeP 




Figure 1: Band structure of LaOFeAs (left) and LaOFeP 
(right). Red and green bands have mainly Fe-3d and 
Pnictogen/Oxygen-p character respectively. The colors are 
guides to the eyes in order to facilitate comparison of the rel- 
ative positions and widths of the bands ( see also blue arrows). 



III. LOW-ENERGY BANDSTRUCTURE AND 
FERMI SURFACE: GREAT SENSITIVITY TO 
FE-PN DISTANCE 

In this section, we focus on the very low-energy 
band structure and Fermi surface properties of the 
LaOFeAs and LaOFeP compounds, as depicted on Fig. [2] 
and Fig. [3] (calculated for the experimental crystal struc- 
ture as given in Table [J). 

In both compounds, the Fermi level is crossed by five 
bands with mainly Fe character. Both materials present 
two electron pockets centered at M and two hole pock- 
ets centered at T with mainly d xz and d xz / yz character. 
(Our choice of the local coordinate system centered at 
Fe is such that the x/y axes point towards the nearest- 
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Figure 2: Low energy region of LaOFePn' s band structures 
with Pn= As (up panel) and P (lower panel). 



neighbour Fe sites. With this choice, the d xy orbital 
is the one that points to the Pn atom while d x 2_ y 2 
points towards the nearest-neighbor Fe.) The signifi- 
cant difference between these compounds comes from the 
third pocket centered along the T — > Z direction. In 
LaOFeAs this third pocket has mainly d xy character and 
presents almost no dispersion along the z-direction (it is 
therefore two-dimensional (2D) in nature). In contrast, 
in LaOFeP it consists of a dispersing 3D pocket with 
mainly d 3z 2_ r 2 character. 

This difference in the Fermi surfaces may explain the 
fact that LaOFeAs is magnetically ordered with better 
nesting properties than LaOFeP , for which no magnetic 
ordering has been reported. When the situation with 
the superconductivity of undoped LaOFeP gets eventu- 
ally clarified, this difference in the Fermi surfaces may 
actually prove to be significant also in this respect. 

The explanation for this difference in the low-energy 
band-structure, as obtained from first-principles elec- 
tronic structure methods, lies mainly in its great sensi- 
tivity to the Fe-Pn distance and to the vertical distance 
of the Pn to the Fe plane. The reason for this sensitivity 
is that, close to the Fermi level, there is a sizeable hy- 
bridization between Fe and Pn states as will be shown 
later in details. 

We checked the effect of doping by using the virtual 
crystal approximation. We observe that for the experi- 
mental structures there are no qualitative changes in the 
low energy band-structure for LaOFeAs , only a shift of 
the chemical potential with still three 2D-like pockets 
crossing the Fermi level. In contrast, for LaOFeP , 
the 3D pocket gets filled and almost disappears from 
the Fermi surface in the doped compound. One should 
note that the bandstructure and FS properties obtained 
for the relaxed LaOFeAs structure (in contrast to the 
experimental one) using LDA (also reported by Singh-^ 
) presents features that are very similar to the ones 
described above for the LaOFeP material. This is due 
to the fact that LDA overestimates the Fe-As bonding, 
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LaOFeAs 




Figure 3: Fermi surfaces for LaOFeAs and LaOFeP calculated 
at their experimental crystal structure. 

so that the Fe-As distance (~ 2.4 A experimentally) 
in the relaxed structure becomes as small as in the 
LaOFeP compound, which is around 2.3 A. As a result, 
the LDA relaxed structure of LaOFeAs displays a 3D 
pocket which disappears upon doping. 

In order to understand better the sensitive depen- 
dence of the low-energy bandstructure on the Pn height 
with respect to the Fe plane, we have studied the 
LaOFeAs system by fixing the lattice parameters a and 
c at their experimental values, while varying the vertical 
position of As all the way from the experimental value 
down to a position which is very close to the one that P 
has in the LaOFeP material (see Table H]). 

The results of this study are depicted on FigQ] We see 
from this figure that the 2D pocket with d xy character in 
the 'experimental' LaOFeAs (upper panel) evolves into a 
3D pocket with d 3z 2_ r 2 character as As is moved closer to 
the plane (lowest panel). The circles in the plots indicate 
As contribution (mainly As-p z ) to the d 3z 2_ r 2. 

In a nutshell, one can say that 'LaOFeAs evolves into 
LaOFeP ' as the height of the As to the Fe plane de- 



creases. 
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Figure 4: Sensitivity of the low energy band structure for 
LaOFeAs to d v (As' s vertical distance with respect to the Fe 
plane) for fixed lattice constants a and c. The considered d v 
values are indicated in each plot. 



IV. ANALYSIS IN TERMS OF WANNIER 
FUNCTIONS 

In the previous section we have noted that the low- 
energy band structure, and most notably the topology of 
the Fermi surface depend very sensitively on the internal 
structural parameter z(As). Reducing z(As) amounts to 
decreasing the distance of the As ions to the Fe plane, 
and changing accordingly the angle that the Fe-As bond 
is forming with the Fe-plane. In this section, we ob- 
tain more direct insights into this sensitive dependence by 
constructing Wannier functions, which reveal the strong 
covalency and hybridization effects associated with the 
Fe-As bond. We present calculations done at the experi- 
mental structural parameters. 

We have used two different types of Wannier construc- 
tions. The first one is a maximally localized Wannier 
function construction (MLWF)2i based on an ultra-soft 
pseudopotential approach as implemented in the PWSCF 
codei£. The second one is the N-th Order Muffin Tin 
Orbitals (NMTO) method and downfolding technique^, 
based on an LMTO-ASA approach^. Previous experi- 
ence has shown 2 - that these two approaches give results 
that are quite comparable to each other. 
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1.78 
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2.05 
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1.28 



Table II: Spread of maximally localized Wannier functions, 
both for a model retaining only the Fe-d states and a model 
that also includes As- and O-p orbitals. Note in particular 
the large spread of the xy orbital, which points towards the 
As atom. 



A crucial issue is the choice of the energy range on 
which the downfolding is performed, i.e the set of Bloch 
bands that the Wannier construction aims at reproduc- 
ing. We shall consider two possible choices: a restricted 
energy range covering roughly [—2, +2]eV around the 
Fermi level, corresponding to downfolding on the set of 
2*5 bands with dominant Fe character. We shall call 
this the 'd-downfolding'. We have also consider down- 
folding on a much larger energy range (approximately 
from — 6eV to +2eV) encompassing all bands with Fe-d, 
O-p and As-p character (2 * 11 bands), which we shall 
denote by 'dpp-downfolding'. 

In Table II, we display our results for the spatial ex- 
tension of the various Wannier functions of LaOFeAs ob- 
tained by the MLWF construction, for both choices of 
downfolding, as measured by the spread: Q = \J\ < r 2 > 
— < r > 2 ). The most striking feature is the large val- 
ues of the spreads, in particular of the d xy orbital which 
points towards the As-atoms, when downfolding is per- 
formed on the restricted set of Fe bands. Note that the 
spread of the d xy orbital is actually comparable to the Fe- 
As distance. The spreads are considerably reduced when 
performing the dpp-downfolding involving the larger set 
of bands, as expected. 

In order to visualize this effect, we plot in Fig. UK a) the 
isosurface of the xy-like Wannier function. The 'leakage' 
of the Fe state on the neighboring As atoms is clearly 
visible. In Fig. [UJb) we plot the same isosurface for the 
same d xy orbital, but in the 'dpp' Wannier construction. 
As expected, the orbitals are considerably more localised 
in this case, and the leakage onto the neighboring As 
atoms is suppressed. 

In passing, we note that these observations have direct 
implications for many-body calculations of these mate- 
rials performed within the dynamical mean-field theory 




Figure 5: (a) Fe-xy Wannier function for LaFeOAs, from the 
maximally localized Wannier construction performed for the 
restricted set of Fe-like bands (d-downfolding). (b) Same or- 
bital character, but within the dpp downfolding, retaining the 
full set of Fe-d As-p and O-p bands. 



(DMFT) framework. Because of the spatially extended 
character of the Fe- Wannier functions in the restricted 
d-downfolding, it seems that including only local matrix 
elements of the Coulomb interaction and local compo- 
nents of the self-energy would be a rather poor approxi- 
mation for this 'd-only' model. A full model using dpp- 
downfolding (hence with more local Wannier functions) 
seems like a better starting point for a DMFT calcula- 
tion. This observation may help understanding some of 
the differences recently reported in DMFT calculations 
of LaOFeAa 2 -^. However, another important issue is 
the value of the on-site Coulomb interaction and the effi- 
ciency of screening, which determines whether these ma- 
terials should be viewed as in the intermediate or strong 
correlation regime. 

Finally, in Fig. [6l we display cuts through the isosur- 
faces of some of the Wannier functions obtained within 
the NMTO construction within d-downfolding. More 
specifically, we have chosen to visualise the most ex- 
tended orbital (xy), the least extended one (x 2 — y 2 ), 
and the 3z 2 — r 2 which due to its orientation along the 
z-axis plays a special role. Also, as shown in the previous 



section, the roles of the xy and the 3z" 



orbitals some- 



what interchange concerning the Fermi surface properties 
when the z parameter is varied. 

We display three specific cuts: the first two are parallel 
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to the Fe xy-plane and correspond to (i) the Fe-plane, 
(ii) a plane containing the lower layer of pnictide atoms. 
The last one is (iii) a cut along a plane parallel to the 
z-axis and containing a nearest-neighbour Fe-As bond. 
It thus cuts through half of the La, Fe, and As atoms 
in the unit cell. While the first cut mainly gives the 
orientation of the Fe orbital, the other two are measures 
of the leakage of the iron orbitals onto the As states. The 
cuts are represented in the three columns of Fig. [6j the 
rows show three different Fe-d orbitals: xy, 3z 2 — r 2 and 

2 2 

x z - y*. 

The first panel in Fig. [6] shows the xy-orbital - the one 
pointing towards the As atoms - of the Fe atom at the 
center of the plot. Also clearly visible are the contribu- 
tions of this orbital on the neighboring four Fe atoms, as 
well as parts of the lobes of the four next nearest neigh- 
bor Fe atoms. Note, that the two As atoms do not lie 
in the plane represented here, but 1.32 Aabove and be- 
low respectively. Still, they mediate the hopping to the 
next nearest Fe atoms, which display more extended xy- 
character than the nearest neighbours, consistent with 
the observation of strong nearest neighbor hopping. The 
second plot in the first row of Fig. [6] shows this same Fe 
orbital but represented by an isosurface cut containing 
the As-plane. The contributions of the Fe orbital leaking 
into this plane and in particular onto the two As atoms 
(eye-like features) are clearly visible. Finally, the cut of 
this orbital in the plane containing the z-axis confirms 
the very extended nature of this orbital. 

The second row in Fig. [6] shows the same cuts, but for 
the 3z 2 — r 2 orbital, the third row for the x 2 — y 2 orbital. 
Particularly interesting are the second plots in these rows 
which visualise the hybridisation with the p x ±ip y orbitals 
of the As atoms, respectively. For LaFeOP (not shown) 
a similar picture emerges. 

The overall comparison of the different orbitals con- 
firms the result of Table 1 that the xy orbital is by far 
the most extended one (followed in fact by the yz and xz 
orbitals - not shown) , displaying a tremendous contribu- 
tion on the neighboring As sites, in a pancake-like flat 
shape. Fig. [7] shows the partial densities of states (DOS) 
of the different Fe-d orbitals, in comparison with the As 
partial DOS. The pronounced peaks of the xy and the 
xz/yz densities of states in the region of the As-p bands 
demonstrate once more that the large spatial extension 
of these orbitals is due to their hybridising with the As-p 
states. 

The Wannier analysis can thus be summarized saying 
that it shows the large extension of the d-orbitals when 
a d-downfolding is used, with the xy- [3z 2 — r 2 ] being the 
most [least] extended one. This leakage is consistent with 
the high sensitivity of the low energy electronic structure 
with respect to the z parameter observed in the previous 
section. 



V. WIDTH OF THE IRON BANDS: 
MATERIALS DEPENDENCE 

Having documented in the previous sections (especially 
from a comparison of LaOFeP and LaOFeAs ) the im- 
portance of covalency and of the Fe-Pn distance, we fi- 
nally turn to the influence of structural properties on the 
bandwidth of the iron bands, in a broader context. This 
question has very recently attracted attention, with the 
observation by Zhao et aZ.|27] of an apparent systematic 
correlation between the superconducting critical temper- 
ature of several different As-based pnictides and the angle 
formed by the Fe-As bond with respect to the Fe-basal 
plane. Specifically, these authors observed that the var- 
ious materials under consideration have comparable Fe- 
As distance (bond- length) dFe-As while the in-plane lat- 
tice parameter a decreases as T c increases. Furthermore, 
these authors suggested that this may imply a smaller 
bandwidth W for materials with highest T c , and hence 
that superconductivity is enhanced by increasing corre- 
lation strength (increasing U/W ratio). Note, however, 
that this would mean a smaller bandwidth with decreas- 
ing Fe-Fe distance. 

Denoting by the angle formed by the Fe-As bond 
with respect to the basal plane as sketched in Fig. [5] (so 
that the opening angle at the top of the Fe4 As pyramid is 
#3 = 7T — 29 in the notations of Ref . l27h , the interatomic 
distances characterizing the Fe-As unit and the vertical 
distance of As with respect to the basal Fe plane are given 
by: 




Figure 8: Schematic Fe4As pyramid for LaOFeAs. The 9 and 
83 angles are indicated. 

In order to investigate the dependence of the band- 
width on these structural parameters, we performed the 
following studies, which can be viewed as 'numerical ex- 
periments': 

• (i) Keeping the lattice constants a and c at their ex- 
perimental values in undoped LaOFeAs , we stud- 
ied how the band-structure evolves as the As is 
brought into the plane, starting from its experi- 
mental position down to a value comparable to the 
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Figure 6: Fe orbitals for the 5-band model for LaFeOAs. First row: xy orbital. Second row: Sz 2 — r 2 orbital. Third row: x 2 — y 2 
orbital. In the Fe-plane (first column), in the lower As-plane (second column) and in a plane parallel to the z-axis containing 
one La, O, Fe, As per unit cell (third column). In the first two columns the x-axis is along the diagonal of the plots. Note that 
the x 2 — y 2 orbital has nodes in the latter plane and is thus not shown in this representation. 



height of P in the LaOFeP compound, i.e reducing 
d v or 2 (As) (hence decreasing 9 and increasing 6*3). 
The resulting changes in the low-energy electronic 
structure corresponding to these calculations have 
been reported above in Sec. Hill 

• (ii) Keeping constant the Fe-As distance e?Fe-As = 
2.407 A(as well as c at the same value than above), 



we reduce the lattice parameter a. In so doing, 
the vertical distance d v separating As from the Fe 
plane increases, resulting in a larger angle 6 with 
the basal plane and a smaller angle at the summit 

03 ■ 

• (hi) Finally, we consider the actual experimental 
structure of the parent compounds LaOFeAs , PrO- 
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LaOFeAs 




Energy (eV) 



Figure 7: Partial density of states of the Fe-d orbitals (solid 
lines), compared to the As density of states (dashed lines, 
same curve in all panels). 



FeAs and SmOFeAs. These three materials have 
almost identical Fe-As distance and correspond to 
decreasing values of a, similarly to (ii), but in con- 
trast to this case they also correspond to decreasing 
values of the lattice constant c, with the ReO unit 
getting closer to the FeAs one. 

The results of these investigations are depicted in Fig. 
in which we display the change in both the full band- 
width W of the Fe bands and the bandwidth of the occu- 
pied Fe states W occ (i.e the distance from the bottom of 
the Fe bands up to Fermi level). Both quantities provide 
useful information, but the latter is usually a better as- 
sesment of the relative importance of correlation effects 
(being related to the typical kinetic energy). 

We note the following trends for the three numerical 
experiments above: 

• (i) As d v is reduced for fixed a, a systematic in- 
crease of the bandwidth (and of the bandwidth of 
the occupied states) is observed. This is indeed 
expected: the direct Fe-Fe hopping is unchanged, 
while hybridisation effects with As are increased 
as d v is decreased, hence increasing the hopping 
through As sites. Thus, in this case, the bandwidth 
decreases as #3 decreases (and 9 increases). 

• (ii) In this case where dFe-As is kept constant, both 
the occupied and full bandwidth show a smaller 
effect as the a parameter decreases. Hence, there is 
a slight dependence of the bandwidth on the angles 
6,63 as compared to i). The reason for this is that 
the increased Fe-Fe hopping (due the decrease of a) 




Figure 9: The change AW in the Fe 3d full bandwidth W 
(solid line) and the bandwidth of its occupied part W occ 
(dashed line) with respect to calculated value for LaOFeAs at 
its experimental structure-. These reference values are 4.25 
eV for W and 1.95 eV for Wocc From the top to bottom 
panels we display the results for the calculations i), ii), iii) 
as explained in the text. Left: AW as a function of dv or a, 
as indicated. Right: AW as a function of S3. The structural 
data for Pr and Sm compounds were taken from Ref. [H and 
I29I . respectively. 



is precisely compensated by the decrease of indirect 
hopping through As (due to the increase of d v ). 

• (iii) Finally, for the real materials which also have 
basically a constant dFe-Ast we observe neverthe- 
less a stronger increase of the bandwidth as a is 
reduced as compared to ii) . The reason for the dif- 
ference with (ii) is that the lattice parameter c si- 
multaneously decreases, and indeed the increase of 
the bandwidth comes also from unoccupied bands 
above Fermi-level, reflecting the influence of the 
Re-0 unit becoming closer. Hence, for the three 
materials investigated, the bandwidth increases as 
03 decreases (and 9 increases), in contrast to (i). 
The occupied bandwidth is even less sensitive to 6*3, 
again because of the compensating effect observed 
for (ii). 

The general conclusion of this investigation is that, as 
expected physically, the bandwidth is controlled both by 
the direct Fe-Fe hopping in the plane and by the indirect 
hopping through As (in view of covalency and hybridiza- 
tion effects). Both the parameters a and 9 (or 63) are 
therefore important. The bandwidth can in principle be- 
have differently (increase or decrease) as the opening an- 
gle 03 is decreased. For the three real materials that we 
have investigated, the bandwidth displays a sizable in- 
crease for decreasing 9 3 however, in contradiction to the 
proposal made in Ref. H3- This increase is basically re- 
lated to the decrease of the c-parameter and the increase 
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Fe-Fe hopping (at constant c the bandwidth experiments 
a more slight effect). 

VI. CONCLUSIONS 

The iron oxypnictides family of compounds share a 
similar overall electronic structure. Some important as- 
pects, however, depend very sensitively on small struc- 
tural changes, particularly on changes of interatomic dis- 
tances and bond angles within the iron-pnictide plane. 
Using first-principles full-potential electronic structure 
calculations, we investigate this sensitive dependence, 
contrasting in particular LaOFeAs and LaOFeP , which 
display different transport, magnetic and superconduct- 
ing properties. The width of the Fe-bands is signif- 
icantly larger for LaOFeP , indicating a better metal 
and weaker electronic correlations. These two materi- 
als also present significant differences in their very low- 
energy bandstructure, when calculations are performed 
at their experimental crystal structure. Both materials 
have three hole-pockets around the T point and two elec- 
tron pockets around the M-point. However, one of the 
hole pockets changes from a three-dimensional one to a 
tube with two-dimensional dispersion when going from 
LaOFeP to LaOFeAs. These differences are due to the 
shorter Fe-Fe distance and to the shorter distance of the 
pnictide to the iron-plane in LaOFeP. We have shown 



that the low-energy bandstructure of LaOFeAs evolves 
towards that of LaOFeP as the As atom is lowered from 
its experimental position to a position closer to the Fe 
plane. We argue that the physical origin of this sensi- 
tivity to the iron-pnictide distance is the covalency of 
the iron-pnictide bond, leading to strong hybridization 
effects. To illustrate this, we have constructed Wannier 
functions which are found to have a large spatial exten- 
sion when restricting the energy window to the bands 
with dominant iron character. Finally, we have clari- 
fied how the bandwidth changes as one moves along the 
rare-earth series in ReOFeAs. Since the Fe-As distance 
remains essentially constant, the Fe-Fe distance contracts 
while the vertical distance to the Fe plane increases, re- 
sulting in a compensating effect and rather small changes 
of the bandwidth. A slight increase of the bandwidth of 
unoccupied states is observed, which is actually associ- 
ated with the decreasing distance between the Fe-As and 
Re-0 planes. 
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